Revealing the physics of r-modes in low-mass X-ray binaries 
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We consider the astrophysical constraints on the gravitational-wave driven r-mode instability in 
accreting neutron stars in low-mass X-ray binaries. We use recent results on superfluid and super- 
conducting properties to infer the core temperature in these neutron stars and show the diversity of 
the observed population. Simple theoretical models indicate that many of these systems reside inside 
the r-mode instability region. However, this is in clear disagreement with expectations, especially 
for the systems containing the most rapidly rotating neutron stars. The inconsistency highlights the 
need to re-evaluate our understanding of the many areas of physics relevant to the r-mode insta- 
bility. We summarize the current status of our understanding, and we discuss directions for future 
research which could resolve this dilemma. 
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What limits the spin rate of a neutron star? Given 
that the fastest rotating neutron star (NS) in a low-mass 
X-ray binary (LMXB; binary system in which X-rays are 
produced when matter is accreted onto the NS from a 
low-mass stellar companion) and the fastest radio pulsar 
have spin rates which are significantly below the centrifu- 
gal break-up limit, it is natural to ask whether a physical 
mechanism prevents further spin-up during the evolution 
of these systems. The issue is challenging because of com- 
plex, often poorly understood, physics. One possibility is 
that the emission of gravitational radiation from the NS 
plays a significant role. Alternatively, the answer could 
be due to the detailed nature of the accretion of matter 
(and angular momentum) onto the NS surface. Despite 
its obvious importance, the question remains unresolved, 
with only the current gravitational wave (GW) searches 
and X-ray observations serving as constraints [l|. 

This Letter concerns one of the main mechanisms that 
is expected to affect the spin evolution of an accreting 
star: the instability associated with the r-modes, which 
are a class of oscillations in a star whose restoring force 
is the Coriolis force. The emission of gravitational waves 
can excite r-modes in the NS core and cause the ampli- 
tude of the oscillations to grow. The notion that this 
instability can provide a spin-limit for NSs in LMXBs 
was first discussed in [2| . The r-mode instability is inter- 
esting for many reasons, mainly because the associated 
gravitational wave signal may be detectable with ground- 
based instruments, but also because its understanding 
requires knowledge from a wide range of physics. The 
primary agents that enter the r-mode discussion are (1) 
damping mechanisms related to the standard shear and 
bulk viscosities and exotica like hyperons, quarks, and su- 
perfluid vortices, and (2) the fluid dynamics associated 
with the mode, e.g., nonlinear coupling and saturation 
[3|-l5|. The instability depends primarily on the NS spin 
rate Vs and core temperature T . This leads to an insta- 



bility "window," determined by a critical curve (defined 
by the balance of evolution timescales tqw = Tdamp) in 
the Vs-T plane, inside which the instability is active. So 
far, most studies of the unstable r-modes focused on par- 
ticular damping mechanisms, in order to determine the 
extent to which they can kill the instability or are rela- 
tively unimportant. In several cases, e.g., hyperon bulk 
viscosity, the answer changed as our understanding im- 
proved ja, lD]. For accreting NSs in LMXBs, the general 
view is that the main damping mechanism is related to a 
viscous boundary layer at the crust-core interface [8|, |9[ . 
What has not been appreciated is that this model leaves 
the majority of the observed LMXBs significantly inside 
the instability window: rapidly rotating NSs should not 
possess spin rates at their observed levels. This is the 
primary message of our Letter and one that requires at- 
tention and resolution. 

Previous discussions considered the general region in 
Vs-T where the LMXB population resides [3|, l8| . Here we 
provide detailed estimates of the likely core temperatures 
for each LMXB, accounting for nucleon superfluidity and 
superconductivity at the level indicated by recent results 
from the NS in the Cassiopeia A supernova remnant (the 
youngest NS in the Galaxy) [10|, [ill . This allows us to 
identify specific LMXBs that are most likely to exhibit 
signatures of the r-mode (in)stability. We bring together 
the theoretical models that have been examined in the 
past (involving, e.g., elasticity, exotic states of matter, 
and superfluidity) and demonstrate that, based on our 
current state of knowledge, they all fail to explain the 
observed systems. This dilemma is irrespective of the 
Cassiopeia A superfluid results. We discuss the uncer- 
tainties and outline where advancements can be made. 

Neutron star core temperatures.— Let us assume that 
the r-mode instability is active in LMXBs at the level re- 
quired to balance the accretion torque, while the associ- 
ated heating is balanced by neutrino cooling. The accre- 



tion luminosity Lace and NS spin frequency Vg (~ f^s/27r) 
are measured from LMXB observations. Since the NS 
(with mass Af and radius R) is taken to be in spin- 
equilibrium, the spin-up torque from accretion is equal 
to the spin-down torque from gravitational radiation, 
i.e., iVacc = ^GW- We take iVacc = Lacc/^K, where 
f^K = {GM/R^y/^ is the Kepler rotation frequency, and 
use the model of [3|, ll^l to obtain A^gw for a r-mode with 
amplitude a and timescales t for relevant processes. Con- 
sidering a 1.4Msun, 12.5 km NS, the balance yields an 
"equilibrium" r-mode amplitude 
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Our choice of M and R are in line with those used in 



previous r-mode work [12| , as well as those given by the 
Akmal-Pandharipande-Ravenhall equation of state (see 
below). In steady-state, the heat dissipated by damping 
of the r-mode is equal to the energy gain from GW emis- 
sion, i.e., Lhoat = -Lgw, where Lqw = -^gw^s/3, so 
that [H 



-'heat 



iacc^s/Sr^K = 0.065(l^s/300 Hz)La 
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Taking the heat from r-mode dissipation to be lost by 
neutrino emission [Lheat = L^{T)], the core temperature 
T can be inferred. Note that cooling via neutrino emis- 
sion dominates over photon emission at our considered 
temperatures. It is traditional to assume that the NS 
cools by the modified Urea neutrino emission process, 
which has a luminosity 
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L^^ « 7.4 X 10^1 ergs s"! (T/IO^ K)^ 



(3) 



Setting Lheat equal to L^^ yields the core temperature. 
Figure [1] shows Lhcat and L^^ , with their intersection 
indicating the core temperature for each LMXB. It is 
worth noting that the heat associated with the unstable 
r-mode [Eq. ^] corresponds to ^ 10 MeV per accreted 
nucleon, compared to the ~ 1 MeV from nuclear burning 
in the deep crust [l5[ . Even if we assume nuclear heating 
instead of an unstable r-mode, the effect on the inferred 
T is only at the < 30% level because of the strong tem- 
perature scaling in Eq. ([3]). 

The above estimates assume normal nucleons in the 
stellar interior. It is expected that neutrons are super- 
fluid and protons are superconducting in the NS core [16| . 
The measurement of rapid cooling of the Cassiopeia A NS 
111 Il7| gives the first direct evidence for the existence of 
superfluid components and constrains the critical tem- 
peratures for the superfluid transition Ten and Tcp, i.e., 
Tcn,max « (5 - 9) X 10^ K and Tep - (2 - 3) X 109 K 



lot lllj . Superfluidity has two important effects on neu- 
trino emission and cooling: (1) suppression of emission 
mechanisms, like the modified Urea process, that involve 
superfluid constituents and (2) enhanced emission near 
the critical temperatures due to Cooper pair formation 
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FIG. 1: Heat generated by damping of r-modes Lheat com- 
pared to the neutrino cooling luminosity L^ as a function of 
NS core temperature T. The thin horizontal lines are Lheat 
for known LMXBs computed using their flux, distance, and 
spin frequency from [2g| and Eq. ([2]). The long-dashed line is 
the modified Urea luminosity Lj^''^. The triangles and starred- 
triangles indicate the intersection of Lhcat and L^^ , which de- 
termines T for each LMXB and short recurrence time LMXB 
(sLMXB). The thick solid lines are L^^ with Tcn^max = 5.6 
and 9 x 10* K, and the squares and diamonds are the inferred 
T (from I/hcat = Lf/) for each source. The short-dashed and 
dotted lines are approximate fits to Lf,^ in the strongly super- 
fiuid and in the non-superfluid neutron regimes, respectively. 



[18| . We use the results of [l9| to calculate the neutrino 
emissivities due to the modified Urea process, accounting 
for superfluid suppression, and the Cooper pair forma- 
tion process. We take Tcp = 2 x 10^ K and Tcn(p) to be 
approximately given by model (a) of [ll|. The neutrino 
luminosity Lf,^ is then obtained by integrating the emis- 
sivities using a stellar model based on the APR EOS 
with M = 1.4Ms„n and i? = 12 km ^. The results 



1.4Msun and i? = 12 km 
presented here do not depend strongly on the assumed 
stellar mass 1G|.Ii1|. 

Figure [2 shows T^f ^jth Tc„,,nax = 5.6 and 9 x 10^ K. 



At T > Ten, max, the supprcssion of the modified Urea 
process by the superconducting protons yields Lf,^ < 
L^^ . Since cooling is less efficient, the inferred core tem- 
peratures are higher than those obtained from Eq. ^. 
At T < Ten, max, neutrino emission is enhanced due to 
Cooper pair formation, and the cooling is more efficient, 
which results in a lower inferred T. It is noteworthy 
that, for Ten, max ^ 8 X 10^ K, a unique T does not 

heat- For example, for 



exist for a range of Lf,^ = L 
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en, max 



5.6 X 10 K, there can be a factor of two differ- 
ence in the inferred T when the observed accretion lumi- 



nosity iacc ~ (3 - 10) x lO^^^ ergs s-i(600 Rz/vs). Inter- 
estingly, there are five LMXBs that show short recurrence 
times between multiple X-ray bursts due to nuclear burn- 
ing of accreted matter [20[ . These sources have accretion 
luminosities within this range, and thus their higher tem- 
peratures could perhaps be responsible for their distinct 
bursting behavior. We also note that there is a branch 
of Lf,^ that could produce LMXBs which increase in lu- 
minosity even though their temperatures are decreasing. 
Finally, we find that, in the temperature regime (T <^ 
Ten, max) whcre both protons and neutrons are strongly 
supcrfluid, the neutrino luminosity is [ISj 
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while in the temperature regime (Tcn,max ^ T' *C Tcp) 
where protons are superfluid and neutrons arc normal, 
the neutrino luminosity is 



4 X lO'^^ ergs s-i[log(r/10* K) 



21 



(5) 



Figure [T] also shows L"p^^ and Lf,^^ . By setting Lhcat 
equal to L"P^^ or Lfp^ , we obtain core temperatures 
which approximate the ones illustrated in Fig. [2] 

Physics of the instability window.— Figure [5] shows the 
core temperature (inferred from either Lheat = L^^ or 
-^hcat = Lf^) and spin frequency for each LMXB. Since 
superfluidity suppresses damping mechanisms like hy- 
peron bulk viscosity and alternative mechanisms like mu- 
tual friction are too weak (see below), the consensus view 
is that the viscous boundary layer at the crust-core in- 
terface is the primary damping agent. It is clear that 
a large number of LMXBs are in the unstable region 
(above the rsv-curve) unless the damping is described 
by a rigid crust model (rBL-curve) [8|]. However, a rigid 
crust is completely at odds with expectations. In the 
fast systems, the Coriolis force that drives the r-modes 
should dominate the elastic restoring force (/x/rig ^ 10""*, 
where ^ is the shear modulus). "Slippage" between the 
crust and core reduces the damping by a factor >100 
(see Fig. [5]) [9[ . It is also worth noting that the magnetic 
fields in these systems {^ 10* G) are too weak to alter 
the nature of the boundary layer (this requires core fields 
> 10^^ G [2l|). We consider the imphcations of the data 
in Fig. [2] in light of these arguments. 

First, let us assume that the r-modes are unstable. 
One might expect the unstable systems to exhibit a dis- 
tinctive behavior. An example may be the short recur- 
rence time LMXBs, which would make them interest- 
ing targets for gravitational wave searches; we estimate 
that dissipation from an unstable r-mode can power the 
observed quiescent luminosity of these higher tempera- 
ture LMXBs (c.f. [l3[)- Conversely, the low temperature 
LMXBs may be r-mode stable; this idea is supported by 
the LMXBs SAX J1808.4-3658 and IGR J00291-f5934, 
which have measured spin evolutions that are consistent 
with magnetic dipole losses without gravitational radia- 
tion 22|; note that the low temperature LMXBs could 
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FIG. 2: Neutron star spin frequency i^s and core tempera- 
ture T. The long-dashed-dotted lines are the shear viscosity 
instability curve (where tqw ~ tsv), the short-dashed-dotted 
lines are the rigid boundary layer instability curve (where 
T"Gw = Tbl), and the solid lines are a more realistic, elas- 
tic boundary layer curve (where tgw = 10^ tbl). Top panel: 
The triangles are LMXBs whose T are derived from their ob- 
served I/acc and assuming L^^ for cooling. Bottom panel: 
The squares are LMXBs and diamonds are short recurrence 
time LMXBs (sLMXB) whose T are derived from their ob- 
served I/acc and assuming Lf,^ for cooling. 



have even lower temperatures, if, e.g., fast neutrino cool- 
ing processes operate in these sources [l8|. Consider a 
NS that enters the unstable region. The r-mode then 
grows rapidly to an amplitude such that nonlinear cou- 
pling to other modes causes the instability to saturate Q ; 
the saturation amplitude is expected to be much larger 
than that required for spin-balance [c.f. Eq. ([T])]. The 
subsequent evolution is likely to be quite complex [5|. 
In principle, the NS will heat up and spin-down, and 
the LMXB should leave the instability window in a time 
much shorter than the age of the system [23| . Therefore 
the observed LMXBs should all be stable, which contra- 
dicts the data in Fig. [51 Most importantly, all reasonable 
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FIG. 3: Three scenarios that could explain r-mode stability 
in the observed LMXBs. Left panel: Crust mode resonance 
at 600 Hz. Middle panel: Superfluid hyperons (based on [3 
with X ~ 0.1). Right panel: Strong vortex mutual friction 
(based on the strong/weak superfluidity models from [23| with 
B ~ 0.01). The dashed lines indicate the break-up limit. 



evolutionary scenarios [5|, |23| predict maximum NS spin 
rates that are far below those observed. 

For r-mode stability, a revision of our understanding 
of the relevant damping mechanisms is required. We 
consider possible resolutions, starting with the viscous 
boundary layer. The crust-core transition may be more 
complex than has been assumed thus far. This should 
be expected given the presence of a type-II supercon- 
ductor in the outer core of the star [l6| . The details of 
the transition are likely to strongly affect the instability 
window, but the problem has not attracted real atten- 
tion. Crust physics may also be vital. There may be 
resonances between the r-mode and torsional oscillations 
of the elastic crust [9|. Such resonances would have a 
sizeable effect on the slippage factor, leading to a com- 
plicated instability window. Figure [3] gives an example; 
the illustrated instability window has a relatively broad 
resonance at 600 Hz, which is the typical frequency of the 
first overtone of pure crustal modes. Although our ex- 
ample is phenomenological (c.f. [Sj), it suggests that this 
mechanism may explain the stability of LMXBs. Realis- 
tic crust models are needed to establish to what extent 
this is viable. 

Another possibility is an instability window that in- 
creases with temperature [2J]. If this is the case, then 
LMXBs may evolve to a quasi-equilibrium where the r- 
modc instability is balanced (on average) by accretion 
and r-mode heating is balanced by cooling (as in our tem- 
perature estimates). This solution is interesting because 
it predicts persistent (low-level) gravitational radiation. 
Figure [3] shows a model using hyperon bulk viscosity sup- 
pressed by superfluidity. However, this explanation has 
a major problem. We must be able to explain how the 
observed millisecond radio pulsars emerge from the ac- 
creting systems. Once the accretion phase ends, the NS 
will cool, enter the instability window, and spin down to 
^ 300 Hz (see Fig. |3|). In other words, it would be very 
difficult to explain the formation of a 716 Hz pulsar j25| . 

A more promising possibility involves mutual friction 



due to vortices in a rotating superfluid. The standard 
mechanism (electrons scattered off of magnetized vor- 
tices) is too weak to affect the instability window [26|. 
However, if we increase (arbitrarily) the strength of this 
mechanism by a factor ~ 25, then mutual friction dom- 
inates the damping (see Fig. [3|). Moreover, this would 
set a spin-threshold for instability similar to the highest 
observed Vs and would allow systems to remain rapidly 
rotating after accretion shuts off. Enhanced friction may 
result from the interaction between vortices and proton 
fiuxtubes in the outer core, as proposed in a model for 
pulsar free precession [27| . This mechanism has not been 
considered in the context of neutron star oscillations and 
instabilities, but it seems clear that such work is needed. 

In summary, we considered astrophysical constraints 
on the r-mode instability provided by the observed 
LMXBs. Having refined our understanding of the likely 
core temperatures in these systems using recent supcr- 
fiuid data, we showed that several systems lie well in- 
side the expected instability region. This highlights our 
lack of understanding of the physics of the instability and 
the associated evolution scenarios and at the same time 
points to several interesting directions for future work. 
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